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We show that the recently developed optical lattices with Peierls substitution - which can be modeled as a 
lattice with a complex tunneling coefficient - may be used to induce controllable quantum transport of ultra-cold 
atoms. In particular, we show that by ramping up the phase of the complex tunneling coefficient in a spatially 
uniform fashion, a finite quasi steady-state current (QSSC) ensues from the exact dynamics of non-interacting 
fermions. The direction and magnitude of the current can be controlled by the overall phase difference but not 
the details of the ramp. The entanglement entropy does not increase when the QSSC lasts. Due to different 
spin statistics, condensed non-interacting bosons do not support a finite QSSC under the same setup. We also 
find that an approximate form of the QSSC survives when perturbative effects from interactions, weak harmonic 
background traps, and finite-temperature are present, which suggests that our findings should be observable with 
available experimental capabilities. 
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There has been growing interest in studying quantum trans- 
port phenomena of ultra-cold atoms. For instance, recent ex- 
periments HJISl have demonstrated how to induce a mass cur- 
rent using different confining potentials and drives, while the- 
oretical studies [6-10] have explored interesting phenomena 
that may soon be tested experimentally. A mass current has 
been usually induced by a displacement or distortion of the 
confining potential 0}|4], a quench or imbalance of the den- 
sity or interactions within the system [6-8 , 10], or connecting 
to external reservoirs ]5]|9). At first glance it may then look 
counter-intuitive that one may induce a current in an isolated 
atomic cloud without either distorting the trap potential or by 
introducing any imbalance in its parameters. Here, instead 
we show the possibility of driving a current under those con- 
straints by using the recently developed techniques of optical 
lattices with artificial (synthetic) gauge fields lfTTlfT2ll . 

Synthetic gauge fields for ultra-cold atoms have suggested 
interesting phenomena such as inducing vortices without ro- 
tating the atomic cloud 1131 , oscillations of atomic clouds 
driven by synthetic electric fields fl4l . and many others. 
Ref. ifTTl shows how a combination of an radio-frequency 
field and a Raman field can generate a lattice potential with 
the feature that atoms acquire a Berry phase when tunnel- 
ing to a different site. This can be modeled as a lattice with 
the Peierls substitution fl5l , where the tunneling coefficient is 
generalized to a complex number ie 1 ^ . Ref. Ifl2l also presents 
optical lattices with complex tunneling coefficient by tempo- 
rally modulating the lattice with designed patterns. As demon- 
strated in Refs. fTT). t and <\> can be tuned separately in exper- 
iments. Here we generalize the concept of Refs. lfTTl[T2l and 
use a micro-canonical formalism Q [16) to study the atomic 
dynamics when is a function of time. We first consider non- 
interacting particles, an important but difficult situation to re- 
alize for electronic solid-state systems which can instead be 
easily accomplished in cold-atom systems flTl . This allows 
us to study the role of interactions in an unambiguous way. 



Then perturbative effects of interactions, finite temperatures, 
and weak harmonic background traps will be studied. 

Our major findings are: i) By ramping up the phase, <j>, as 
a function of time t in a spatially uniform fashion, an atomic 
mass current can be induced in a closed isolated system with 
or without a weak harmonic trap potential, ii) The direction 
and magnitude of the mass current can be controlled by the 
displacement of <j>. Hi) For non-interacting fermions, a quasi- 
steady-state current (QSSC) and its approximate form persist 
for a period of time depending on the system size. A conden- 
sate of non-interacting bosons, in contrast, does not support 
a finite QSSC. iv) The fermionic QSSC exhibits no memory 
effects, which means that the magnitude of the QSSC is insen- 
sitive to the details of the ramp, v) The entanglement entropy 
remains constant when the QSSC lasts, and finally vi) The 
QSSC (and its approximate form) survives when perturbative 
effects from interactions and finite temperatures are consid- 
ered. All these predictions can be tested with minimal modi- 
fications of available experimental techniques lfTTlfT2l . 

Following Ref. fTD . we consider non-interacting ultra-cold 
atoms 1 18 1 loaded into a finite one-dimensional optical lattice, 
which may be modeled by the Hamiltonian 

JV-l 

H = -J2 [t^c]c J+1 + H.c] (1) 

3=1 

with open boundary conditions. Here t and are the ampli- 
tude and the phase of the hopping coefficient and should be 
tunable separately IfTTIl . c] (cj) is the creation (annihilation) 
operator of site j, and N is the lattice size. We assume <j> and 
t are uniform across the whole lattice. The unit of time is 
chosen as to = H/i. 

The system can be driven out of equilibrium by varying <fi 
over time. The current flowing from the left half (j < N/2) 
to the right half (j > N/2) is given by I — -dN L /dt, where 
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Figure 1: Schematics of the setups of fermions for (a) phase-induced 
transport and (b) an effective model for describing (a). Here grey dots 
imply that particles may be in a superposition of quantum states. 

Nl = Efll(4°j)- Explicitly, 

7 = 2Im[ie^(cj v/2 c A72+1 )]. (2) 

The current can be measured by the protocol of Refs. 0171. 
The time evolution of the single-particle correlation matrix 
Cij(t) = (c\(t)cj(t)) follows the equations of motion 
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ie^c h3+1 (t) - ie-Wcij-iit). (3) 

For a closed system as the present one, = if i, j < 1 
or i, j > N. We first consider N p non-interacting fermions 
loaded into the ground state of the lattice with = initially. 

We then let the phase (f>(t) rise according to (/>(£) = 
(t/t^)4> m when t<t$ and </>(t) = (j> m when t > t$. Note that 
4>(t) is spatially uniform across the whole lattice at each time 
t. Figure[T](a) illustrates the setup. The unitary transformation 

= Eti U lk b k with U jk = y^sin (^) e^-W 

leads to if = £f =1 ^fc&fe 6 * with E k = -2t cos (7^1) ■ We 
remark that Ujk depends explicitly on <f> while Ek is indepen- 
dent of <fi. In the energy basis, the initial state of fermions 
corresponds to (b\by) = 5kk'0(N p — k) with cf> = 0, where 
0(x) = 1 if x > and 6(x) = otherwise. The correspond- 
ing Cij (t = 0) can be inferred from the unitary transformation. 

In addition to the current, we also evaluate the entangle- 
ment entropy, s, between the left half and the right half. The 
entanglement entropy for a bipartite system with noninteract- 
ing fermions has been discussed in Ref. |[T9l . Here it can 



be calculated as follows. We define My(t) 



;(*) ^ 



1 < i, j < (N/2) as the sub-matrix of the correlation ma- 
trix for the left half with the eigenvalues Vj (t) . Then 

N/2 

s(t) = -!>(*) log Vj (t) + (1 - Vj (t)) log(l - Vj (t))}.(4) 



For electrons modeled as noninteracting particles flowing 
through a quantum point contact 030, the rate at which s(t) 
increases depends on the transmission coefficient, which mea- 
sures the probability that a particle can tunnel from the left 



half to the right half when the systems exhibit a steady or 
quasi-steady current. 

The current and entanglement entropy of a system with 
N = 512 and N p = 256 are shown in Figure|2]for tf> m = n/2 
and (3/2)7r with t$ — 10to and 20to- We consider a finite ra- 
tio of Np/N as discussed in Ref. [20|. Due to the finite size of 
the system, the wavefunctions will interfere with themselves 
after some revival time, as one can see on Fig. [2] Throughout 
the paper we will focus on the physics before this revival time, 
which we found to be linearly proportional to the system size. 

A key point is that even when the lattice, initial density dis- 
tribution, and (f)(t) are all spatially uniform, there can be a 
current flowing across the middle of the lattice. The mech- 
anism behind this current is because the initial ground state 
for 4> — is not the ground state for cf> > so the spreading 
of the wavefunction leads to a mass current. One can clearly 
see that for the currents there is a plateau for each curve in a 
given time interval. This plateau resembles the quasi-steady- 
state current (QSSC) - studied in Refs. |j7] [16] - induced by 
an inhomogeneous bias or density. We found that for a fixed 
ratio of N p /N, the duration of the QSSC scales linearly with 
N and its magnitude is the same. For a fixed N, the dura- 
tion of the QSSC is the same as N p varies but the magnitude 
decreases in a non-linear fashion as N p decreases. The inset 
of Fig. [2] (a) shows the currents for {N,N p ) = (512,256), 
(512, 128), and (256, 128) with </> m = (3/2)tt and t = 10i 
and demonstrates those features ETl . 

Importantly, the direction and magnitude of the QSSC can 
be controlled by 4> m for fixed N and N p . We also remark that 
the QSSC lasts when d(j>/dt = until the revival time and 
this is partially due to the fact that there is no dissipation in 
an isolated non-interacting system. Another important feature 
is that the magnitude of the QSSC for fixed N and N p is in- 
sensitive to different ways of turning on the phase <fi(t) 11221 . 
Our findings suggest that the system exhibits no memory ef- 
fect when the QSSC has been developed and maintained. 

The entanglement entropy remains the same when the 
QSSC lasts, which is consistent with previous studies with 
different setups [ 19] in the limit where the transmission co- 
efficient is equal to one due to the uniform tunneling coeffi- 
cient. One may understand this behavior by the full count- 
ing statistics which models the particles semi -classically with 
probability T (equal to the transmission coefficient) of tunnel- 
ing through the middle and with probability 1 — T of being 
reflected |fl9l . The entanglement entropy increases because 
this probability distribution introduces correlations between 
the right half (with tunneling particles) and the left-half (with 
reflected particles). For a uniform lattice without any scatter- 
ing potential, however, the transmission coefficient is T = 1 
||23l , so particles tunnel through the middle without building 
any correlations between the two halves of the lattice. 

Let us now analyze the effect of ramping up <fi from dif- 
ferent initial values. The corresponding initial states are de- 
termined by diagonalizing the Hamiltonian ([TJ with different 
(j)(t = 0) and constructing the ground states. The currents 
for different initial and final values of 4>(t) are shown in the 
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Figure 2: (Color online) The fermionic current (a) and entanglement 
entropy (b) from phase-induced transport. We show <f) m = tv/2 in 
black and <j> m = (3/2)n in red. Solid and dashed lines correspond 
to t$ = 10to and 20io, respectively. Here N = 512 and N p = 256. 
The inset of (a) shows how the quasi steady-state current varies as 
N or N p changes. The black solid line, red solid line, and black 
dashed line correspond to (N,N P ) = (512,256), (512, 128), and 
(256, 128) with m = (3/2)w and ^ = 10. The inset of (b) shows 
the currents by ramping up (f> from to 7r (black dashed line), from 
7r/2 to tv (red solid line), and from n to 3n/2 (black dotted line) for 
N — 512 and N p — 256 with t$ = 10to- All three curves coin- 
cide. The current and entanglement entropy of the effective model 
(]} with N = 512, N p = 256, and = 10 are shown in (c) and (d). 
The black solid lines, red solid lines, and black dashed lines corre- 
spond to the three selected Vo(t) shown in the inset of (c). Note that 
Jo* Vo(t)dt are the same for all three cases. 



inset of Fig. |2](b). We deliberately choose (<f)(t = 0),c/> m ) 
as (0, 7r/2), (tt/2, it), and (tt, 3tt/2) with the same ramp time 
= 10to- All three curves coincide. This can be explained 
by the fact that in quantum mechanics, an overall constant 
phase change of the wavefunction does not introduce observ- 
able effects. Therefore the dynamics only depends on the dis- 
placement of the phase. One may notice that the QSSC from 
<fi(t = 0) = to (f) m = 3tt/2 has the opposite sign of that 
from to 7r/2. To observe this reversal of the direction of the 
QSSC, one has to ramp the phase by the full amount. Simply 
going from ir to 3ir /2 does not show this reversed current. The 
ability to reverse the direction of the current is because the fac- 
tor e 1 ^ in the Hamiltonian ([T| selects the preferred direction of 



be reversed. One may also see this reversion of the current 
from the expression of I shown in Eq. Q. This relation also 
implies that there is no current when <fi m = <f>(t = 0) + ir, 
which we have verified numerically. 

To understand why varying <p(t) can drive a QSSC, we 
resort to a classical analogy. Using E = —dA/dt and 
/ = V/R, one expects that the instantaneous current is lin- 
early proportional to d<j>/dt, where cj> = J A ■ dl. Here E 
is an electric field causing the voltage difference V and R is 
the resistance of the classical system. We thus assume that 
d(f>/dt induces a voltage difference V between adjacent sites 
j and j + 1 and model the quantum system by the effective 
Hamiltonian 



H eff - ~t 



N-l 



)+^V j (t)c}c j . (5) 



3=1 



Here Vj(t) = if t < 0, Vj(t) = V (t){N - j) during < 
t < ts, and then Vj(t) = for t > t^. In other words, the 
lattice is tilted during < t < t$ with a potential difference 
Vo(t) between adjacent sites. In the effective model |5]) the 
tunneling coefficient is real and time-independent. Figure [T] 
(b) illustrates this effective model. 

Figure [2] (c) and (d) show the current and entanglement en- 
tropy from the effective Hamiltonian |5| for several functional 
forms of Vo(t) shown in the inset. Two important features im- 
mediately connect the effective model and the phase-induced 
transport: i) a QSSC and ii) a constant entanglement entropy 
when the QSSC persists. Therefore it is the potential differ- 
ence due to d(f>/dt during tj, that drives the QSSC and the 
QSSC can continue after d<f>/dt — or the potential differ- 
ence vanishes. Another important feature is that for differ- 
ent functional forms of Vb (i) with the same Jo* V (t)dt, the 
magnitude of the QSSCs are extremely close to each other. 
This again exhibits the lack of memory effects in the trans- 
port of non-interacting systems because the QSSC is con- 
trolled by the area enclosed by Vo(i), not the detailed form 
of Vo(t). However, our results bear a stark difference between 
the classical and quantum systems: For the quantum system, 
the QSSC is controlled by the total difference of (j>, not its time 
derivative. 

In Ref. it was shown that condensed noninteracting 
bosons cannot support a QSSC when the system is subject 
to a sudden density imbalance. For bosons the Hamilto- 
nian and the equations of motion are identical to those of 
fermions. At zero temperature, all non-interacting bosons 
occupy the lowest energy state so the initial condition is 
(b\.bk>) — N p Skk'3k,i with = 0, where k = 1 denotes the 
lowest energy state. We did not find any plateau in the bosonic 
current, namely initially condensed non-interacting bosons do 
not support a QSSC although a time-dependent <f>(t) can still 
drive a finite current E4l . 

Finally we consider other effects that may arise in ac- 
tual experiments. The first case involves two-component 
fermions with weak onsite repulsive interactions of the form 
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H, 



where n,- 



1,2 for 



the two species, and a is the opposite of a. The equa- 
tions of motion with the standard Hartree-Fock approxima- 
tion have been reported in Ref. (20). We consider an ini- 



tial state of the approximate Hamiltonian for the species 

,-r- U — —TX^"^ 1 I J ^ . . . _L TT „ \ -L TTlfr,-\ m . whoro 



a: H a = 



{n^rijtj, where 

= (1/N) l ( n js)- The unitary transformation Cj CT — 
Sfe Ujkbka that diagonalizes H a can be found and the initial 
state corresponds to l/b\J>k'<r') = 5kk'5 acr '0{N p — k), where 
we assume that there are N p particles for each species. Then 
the system evolves according to the Hamiltonian H+Hjj with 
a linear rising <fi(t) for t < t<f, and <fi = (j> m for t > t^. We 
also assume that (cj^c^) = (cj^Cjg-) for all t. Figure 3 (a) 
shows the currents for a = 1 from the case with U/t — land 
from the non-interacting case with N = 512 and N p = 256. 
One can see that for weak interactions where the mean-field 
(Hartree-Fock) approximation is reasonable, the current is not 
affected qualitatively by the presence of the repulsive interac- 
tions. We remark that the initial state is not the ground state 
of the full mean-field Hamiltonian H + Hjj. Therefore even 
if 4>(t) = for all t, there is still a tiny current (~ 10 _9 iQ ) 
induced by the difference between the Hamiltonians. This is, 
however, a negligible effect. 

Other two possible effects involve the background har- 
monic trap for holding the atomic cloud and finite temper- 
ature. We study them for the non-interacting single-species 
fermions case. The harmonic trap potential introduces a term 
X^(l/2)r7iw 2 d 2 (j - N/2 - l/2) 2 c]c J to the Hamiltonian, 
where m is the mass of atoms, ui is the trap frequency, and 
d is the lattice constant. We let y = (l/2)muj 2 d 2 ft. In the 
presence of a harmonic trap, the density is higher at the mini- 
mum of the trap potential and is lower towards the trap edge. 
If the density reaches rij — 1 at the center (j around N/2), 
that region will become a band insulator and no current can 
flow through it in a one-band model. Therefore N p /N and 
y cannot be too large or the insulating state at the center will 
emerge and block any current. We found that for a lattice with 
N = 512 and N p — 128 held by a harmonic trap potential 
with y = 1 x 1CP 4 , the maximal density at the center is less 
than one. Fig.[3](b) shows the density profiles with and with- 
out a harmonic trap potential. The QSSC can still be observed 
at T = when a weak harmonic trap potential is present, 
albeit it lasts for a shorter time period compared to the case 
without the trap potential as shown in Fig.[3](c) and (d). 

To study finite-temperature effects, we first diagonalize 
the <fi — Hamiltonian with or without the trapping poten- 



tial to obtain H to t = Ylk=i -^'ifc(^ t )L^fe- The initial condi- 
tion is chosen as ((&*)£&£/) = f(Ef.)Skk', where f(x) = 
(exp[(x — /i)/fcsT] + 1) _1 is the Fermi distribution function. 
The chemical potential is determined by N p — Ylk=i f(^k)- 
One can use the unitary transformation that diagonalizes the 
Hamiltonian to obtain Cjj(i = 0) and follow the equations of 
motion as <f>(t) turns on. The currents from the initial states 
at T = At/ ks with and without a harmonic trap potential are 
shown in Fig. [3](c) and (d) for 4> m — tt/2 and 37r/2 with 
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Figure 3: (Color online) The phase-induced currents under differ- 
ent conditions: (a) Effects of interactions. We show the currents for 
one species of a two-component Fermi gas with onsite interactions 
U/t — 1 (solid lines) and without interactions (dashed lines). Here 
4> m = ty/2 (black) and m = 3tt/2 (red) with t$ = Wt , N = 512, 
and N p = 256. (b) The initial ground-state density profiles for dif- 
ferent confining potentials and temperatures: (1) and (2) correspond 
to the system without a harmonic potential, while (3) and (4) corre- 
spond to the system with a harmonic potential with y = 10~ 4 . The 
temperatures of (1) and (3) are T — while those for (2) and (4) 
are T = 4f/fc_B. (c) and (d) show the currents from the four cases 
shown in (b). The maximal phase in (c) is <f> m = tv/2 and that in 
(d) is 3n/2. The ramp time is t$ = 10to for both (c) and (d). Here 
TV = 512 and N p = 128 for (b), (c), and (d). 



<0 = I0t . One can see that the phase-induced current at 
finite T is smaller when compared to the current at T = 0. 
This implies that the phase-induced current is a phenomenon 
that requires coherence of the initial state. At finite T the ini- 
tial state is mixed according to the Fermi-Dirac distribution 
and the tuning of the tunneling coefficient could not drive a 
current efficiently. Moreover, in the presence of a harmonic 
trap potential at finite T, we observe a slightly tilted regime 
in the current instead of a plateau that defines a QSSC. There- 
fore in experiments with a harmonic trap, there should only 
be an approximate QSSC where the current changes relatively 
slowly as time evolves and this approximate QSSC represents 
the QSSC that is more prominent as T — > 0. Nevertheless, for 
zero or finite temperatures with or without a harmonic poten- 
tial, the direction and magnitude of the phase-induced current 
could still be controlled by the maximal phase <fi m even when 
4>(t) and t are uniform across the whole lattice. 

We have thus shown that controllable phase-induced trans- 
port should be feasible in optical lattices with Peierls substi- 
tutions in realistic experimental conditions. Combined with 
reservoirs [5|, density- or interaction-imbalance J7] |8], and 
other types of proposed devices ll25ll . phase-induced transport 
may find applications in atomtronics 11251 or quantum quench 
dynamics ll26l and provides an interesting example of dynam- 
ical quantum phenomena. 
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